The enormous outer Galaxy H II region CTB 102 
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iHj; ABSTRACT 
> ' 

^ ■ We present new radio recombination line observations of the previously un- 

^ ' studied H II region CTB 102. Line parameters are extracted and physical pa- 

O ' rameters describing the gas are calculated. We estimate the distance to CTB 102 

^ ■ to be 4.3 kpc. Through comparisons with H I and 1.42 GHz radio continuum 

^ ■ data, we estimate the size of CTB 102 to be 100 — 130 pc, making it one of the 

O ■ largest H II regions known, comparable to the W4 complex. A stellar wind blown 

■ bubble model is presented as the best explanation for the observed morphology, 
^ ■ size and velocities. 

Subject headings: ISM: bubbles - ISM: H II Regions - ISM: Individual (CTB 102) 
- Radio lines: ISM 



Introduction 



The radio b right outer Galaxy regio n CTB 102 {i = 93?060, b = +2?810) was first 
cataloged by the IWilson fc BoltonI (Il960l ) radio survey of the Gala ctic plane. The source 
is then mentioned in subsequent Galactic radio surveys including iKallas &: ReichI (jl980l ) 
where it is identified as KR 1. Using radio recombination line (RRL) observations at A3 cm. 
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Lockmaij (119891 ) (H87a, ~3' beam) identified the region as a H II region with a fine brightness 
of Tiine = 18 ± 1.3 mK, a velocity of Vlsr = —61.0 ± 0.8 km s~^ and a full width at half 
maximum (FWHM) of AV = 22.9 ±1.9 km s~^. Radio continuur a images at 1.42 G Hz and 
~ 1' resolution from the Canadian Galactic Plane Survey (CGPS. iTaylor et al.ll2003l ). show 
filamentary structures extending from a bright complex source. From the appearance of the 
structure and a kinematic distance estimate, the region appears to be a very large H II region 
and a major feature in the Perseus spiral arm. 

Yet this major Galactic region is unstudied. Suffering heavy extinction in this direction 
in the Galactic plane, there is no known optical counterpart to CTB 102. The purpose of 
this study is to determine the basic properties of CTB 102, mainly how large in physical 
size it is, and how it influences its Galactic environment. In this paper we present new RRL 
observations towards CTB 102. RRL observations allow direct velocity measurements, and 
along with continuum observations will tell us the density and temperature of any gas in the 
beam at or near thermodynamic equilibrium. 



2. Observations 



RRL observations towards CTB 102 were performed with the 100-m NRAO Green Bank 
Telescope (GBT) during 6 nights in 2006, July 31, August 2-4, 15 & 17. Twelve pointings 
were observed around the CTB 102 complex. These telescope pointings are illustrated in 
Figures [1] and [31 Positions and total integration times for the chosen observations are given 
in Table [1], the observations themselves were divided into 600 s scans. RRLs observed 
were H103a through HI 10a. A 50 MHz bandwidth receiver was used to allow the eight 
recombination lines to be simultaneously observed in the high end of the C-band {u = 
4.8—5.9 GHz). Both polarizations were admitted, and the spectrum consists of 4096 channels 
(0.62 — 0.75 km s~^ per channel). System temperatures ranged from 19 to 26 K, depending 
mainly on the elevation of the source. Average system temperatures for each observation 
are given in Table [TJ As a check of the system's ability to record RRL emission, the bright 
"head" of CTB 102 at i = 93?115, b = +2?835, hereafter CTB102p, was observed for 600 s 
at the beginning of each session. 



3. Data Reduction 



None of the eight 50 MHz bands were seriously affected by radio frequency interference. 
Frequency-switched scans in each linear polarization (YY, XX) were folded individually; 
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since frequency- switching was done in-band (±12.5 MHz), we doubled our effective inte- 
gration time. Each 600 s scan contains 16 spectra. To assure that no hne structure (e.g. 
very extended wings from outflows) is removed, velocities forbidden by Galactic rotation 
(typically -300 km s"^ < Vlsr < -130 km s"^ and +30 km s'^ < Vlsr < +200 km s"^) 
defines a range of baseline velocities. For each of the 16 spectra in every scan, a baseline 
was determined using a fourth-order polynomial fitted to the range of baseline velocities. 
This fitted baseline was then subtracted from every spectrum in each scan. After baseline 
subtraction, the spectra in the individual scans are combined (for every velocity channel, 
intensity values are summed up and then divided by the number of scans) to create averaged 
spectra, one for each line and polarization. At this point in the reduction process, there are 
16 spectra (H103a XX, H103a YY, H104a XX, etc.) for every observation in Table [TJ 

These averaged spectra were regridded to a common channel width (0.67 km s~^) and 
smoothed to a common velocity resolution (1.5015 km s^^). Typically 4-10 of these averaged 
spectra do not show residual wavy baselines in regions of no RRL signal. The only exception 
is HllOo; (polarization XX), which shows a very wide "bump" in the spectrum centered 
around —200 km s~^, extending into the region of the RRL signal. This line and polarization 
is completely excluded from the analysis. To reduce noise, composite spectra, one for each 
filament, are made by combining the averaged spectra that do not show residual wavy 
baselines. The composite spectra typically have a noise level of ATj-ms ^ 1 niK (antenna 
temperature). 



Since RRLs are expected to be quite wide (25-30 km s~^; iLockmanl 119891 ) . a spectral 

resolution of 1.5 km s~^ is unnecessarily fine. A higher S/N can be achieved without loss of 

information by moderate smoothing of the composite spectra, although too much will add 

an artificial width to spectral lines present. We conservatively choose a resolution of 3.0 km 
1 



s 



4. RRL Results and Analysis 

To our final smoothed composite spectra, Gaussians are fit to obtain spectral line pa- 
rameters: hne amplitude (Ti), central velocity (Vlsr) and FWHM (AV). Smoothed spectra 
and the Gaussian fits are shown in Figure [2] and the obtained parameters are presented in 
Table [2J Note that antenna temperature has been divided by the beam efficiency, rji, = 0.92 
for the GET at 5 GHz, to convert to brightness temperature. The uncertainties in Table [2] 
are obtained in Monte Carlo fashion. To the originally obtained Gaussian fit, randomly 
drawn noise from a normal distribution with the same standard deviation as the previously 
obtained AT^^s is added. A new Gaussian is then fitted to the generated spectrum and 
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its parameters stored. After 1000 repetitions, the standard deviation of the distributions of 
spectral hne parameters is combined with the mean fit uncertainties to produce the total 
uncertainties. Figure [3] displays the observed velocities of the filaments/objects, overlayed 
on a CGPS 1.42 GHz continuum image. 



4.1. Velocities 

This study finds that the bright radio source CTB10 2p has a ve l ocity of Vlsr = — 62.66± 



0.05 km s~^, which is 1.7 km s~^ less than the value in iLockmanl (119891 ). The difference is 
probably due to the difference in pointing position, as well as the difference in telescope 
beam size between the studies. The velocity of CTB102p is hereafter referred to as Vref- 
One purpose of this study is to find out the size of CTB 102, i.e. are the filaments and 
objects we see around CTB 102 at the same distance? The velocity gradient (for a flat 
rotation curve with Rq = 8.5 kpc and Vq = 220 km s~^) in this part of the Galaxy is ~ 10 
km s~^ kpc^^. Looking at the last column of Table [21 most filaments and objects have a 
1^ ~ K-e/l ^ 6 km s~^ , so it seems likely that the CTB 102 complex and objects KR 4, 
KR 6 and NRAO 652 are in the same part of the Galaxy. In contrast, WB 43 (—0.8 km s^^) 
is clearly a local H II region, and in no way connected to the CTB 102 complex (Figure |2]). 
Filament 2 is another possible unrelated object (see Sec [7]), its velocity (—78.4 km s~^) 
deviates the most from the velocity of CTB102p, \V — Vref \ ~ 16 km s~^, and is ~ 10 km 
s~^ less than any other velocity observed in this study. Computing the average absolute 
difference for all observations (except WB 43) yields a result of 6.6 km s~^, but excluding 
filament 2 from the average, the value drops to 4.3 km s^^. 



4.2. Line Widths &z Electron Temperatures 



The majority of FWHMs in this study fa ll between 20 km s ^ and 26 km s ^, which is 
typical for dense H II regions (jLockmarull989l ). 



The line width of a RRL depends on mainly two things, thermal motion, which is 
described by the electron temperature (Tg), and microturbulence. Both of these processes 
produce Gaussian profile shapes, so a Doppler temperature (Tp), which provides an upp er 
limit to the true electron temperature, can be adopted as follows (IRohlfs &: WilsonI 120041 ): 



AV = 0.21T; 



1/2 



D 



where is the observed FWHM in km s ^ from Table [2], corrected for the spectral reso- 
lution of 3 km s~^. 



- 5 - 



The LTE electron temperature can be calculated using the line to continuum brightness 



ratio: 



km s 



-1 



6.985 X 10^ 



a(z/, Tg) vGHz 



1.1 



^1.15 



1 



(2) 



where aiv, Tp 



1 and N^^+/N^+ 



i+%e-/%- 

0.08 jRohlfs fc WilsonI bood l The continuum 



brightness temperature Tc is obtained by convolving the CGPS 1.42 GHz mosaic to the 2.5' 
resolution of the GET. The background is estimated individually for each filament (typical 
background estimate uncertainty ~ 0.5 K) and the background is subtracted from the aver- 
age brightness temperature within a 2.5' circle centered on the telescope pointing. Assuming 
temperature follows the power law TjjT^ = {u/u')^^'^ the resulting is scaled using the 
average frequency of the observed RRLs (5.46 GHz). The parameters derived in this way 
for each filament /object are in Table [3l The Tg is quite sensitive to the estimated back- 
ground temperature, especially when the peak brightness temperature is barely above the 
background. Note that To > Tg within the la uncertainty in all cases. To provides a useful 
upper limit to Tg in case the uncertainty in Tg is large. These derived parameters for the 
ionized gas gives an average LTE electron temperature of ~ 6 x 10^ K. This value falls within 
the typical range of electron temperatures for H II regions. 

We note the FWHM of CTB102p determin ed in this study , AV = 18.34 ± 0.08 km 
s~^, is 4.6 km s~^ less than the value for AV in iLockmanl (119891 ). Again, the difference is 
probably due to the difference in beam size and pointing position. 



5. Distance 



With the majority of the filaments' velocities found to be very similar to that of the 
central body, the whole CTB 102 complex subtends about one to two degrees of sky. Its 
distance would be able to tell us how far across its influence extends across the outer disk. 
Here we provide the first ever distance estimate for this H II region. 

To begin with. Figure H] shows a latitude-velocity sl ice in the direction i ^ 93?6, taken 
from the 26-meter telescope survey of Galactic plane H I ( iHiggs fc Tapping||2000l ). Extensive 
intermediate velocity (—90 < Vlsr < —20 km s~^) H I is seen extending from —2° < b < 
+5° in latitude. The H I beyond Vlsr < ~20 km s~^ appears split into three blended 
concentrations. The lower-latitude concentration at 6 = 0° and Vlsr = ~44 km s~^ blended 
together with the concentration near b = +1?4 at Vt,sr = — 68 km s~^ are contemporarily 
taken together to form the Perseus H I arm (see lRobertslll972l . for example). The Outer arm 
is the smaller oblong H I feature up at 6 = +3?6, Vlsr = —82 km s~^, which is also blended 
with the upper-latitude portion of the Perseus arm. 
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As seen in the top panel of Figure HI the contemporary Perseus arm exhibits a very 
substantial tilt in latitude towards negative velocities. This is known as the "rolling" motion 
in this arm. As well, the spl itting of the H I Perseus arm into two concentrations has been 



explained by a spiral shock (iRobertd 119721 ) which is thought to precede the arm. However, 



an alternative explanation for the £, b, V str ucture in the H I suggests that there are actuall y 



three spiral arms beyond —25 km s ^ (e.g. IValled l2008l : iKimeswenger fc Weinbergerlll989l ): 



the Perseus arm (centered at b, Vlsr = 0°, — 44 km s~^), the Cygnus (=Outer) arm 
(+1?4, —68 km s~^), and the "Far" Outer arm up at 6 = +3?6. In this 3-arm interpretation, 
the Perseus and Cygnus arms have more reasonable moderate rolling gradients {dV/db ~ 
—2 km s~^ per degree). The Cygnus arm appears at higher latitudes than Perseus since it 
lies further into the upwardly-warped outer Galactic disk. 

The systemic velocity estimated for CTB 102 in H I (:i; —58 km s^^; see Sec [7]) and 
its latitude {b = +2?8) place it somewhat in the middle of the contemporary Perseus arm, 
in an H I "wall" that connects the two brighter concentrations. In b, V plots (Figure HJ 
top), CTB 102 is situated between the two major H I concentrations, and in the 3-arm 
picture above it is not clear in which spiral arm it resides. In i, b maps (Figure HI bottom) 
a "mushroom-cap" -like feature centered on CTB 102's coordinates {i = 93.2°, b = +2?8, 
Vlsr = ~58 km s~^) appears as a "blister" off of the top of the extensive H I arm at 6 = +2°. 
As well, in an i, V plot of 80° < i < 110° (Figure HI center) this cap is seen as a finger-like 
extension off of the Cygnus H I arm, centered at Vlsr = —57 km s^^ extending down to 
—47 km s^^. It has the appearance of either an upper-latitude extension to the Perseus arm 
(up to b = +3°), or a blended H I concentration (possibly the wall of an expanding bubble) 
extending off the Cygnus arm to more positive velocities. A distance to the Perseus arm in 
this direction then forms a good estimate for at least a lower-limit distance to CTB 102. 

The velocity of CTB 102 {Vlsr — —58 km s^^) indicates a rather large kinematic dis- 
tan ce of 6.8 kpc, assuming a fla t rotation curve with Rq, Qq — 7.6 kpc, 28 km s~^ kpc~^ (e.g. 



see 



Kothes fc Dougherty! 120071 ) . However, this assumes that the systemic Vlsr. derives only 
from the projection onto the line-of-sight of the object's Galactic circular motion. Being near 
a major spiral feature of the Galaxy, CTB 102 is likely influenced by gravitational forces from 
the stellar component of the arm, so its velocity is likely tainted with non-circular motions 
such as those from the "rolling" motions in the arms (velocity gradients perpendicular to 
the plane), and streaming motions from the gravitational influence of density waves. 

A new kinematic-based distance method that accounts for non-circular streaming mo - 



tions due to a two- armed density wave has been developed by [Foster fc MacWilliamd (120061 ). 
The approach is to model the Galactic H I distribution and rotation curve by fitting an em- 
pirical model of Galactic structure and density-wave motions to the observations, rather than 
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just assuming a purely circular model. The model assumes only two arms are present in the 
second quadrant, arranged as a two-armed density wave pattern. We model a direction near 
to 6 = 0° to ens ure that velocity gradients w ith latitude (i.e. rolling motions, which are not 
accounted for in Foster fc MacWilliamsl 2006 ) are not included. The "rolling" gradient of the 



H I Perseus arm observed towards CTB 102 is dV/db ~ —2 km s~^ per degree of latitude; 
corrected for this latitude gradient, CTB 102's systemic velocity is Vlg"^ ~ —52 km s~^. 
The velocity field of the best fitting H I model in the direction i = 92.4°, b = —0.2° is 
shown in Figure O In this direction only the Perseus arm is seen along the line-of-sight. 
The model shows the Perseus Arm gaseous density peak (associated with the spiral shock) 
to be 4.3 kpc distant from the Sun in this direction. The velocity field is ambiguous be- 
tween —55 < Vlsr < —44 km s~^ due to the shock, so the heliocentric distance to CTB 102 
i^LS^R ~ km s~^) is either 4.3 kpc or 5.8 kpc. 

A simple line of reasoning will show that r = 5.8 kpc is probably not the true distance 
to CTB 102. Supposing that it is, then for i = 93° the galactocentric distance to CTB 102 
is -R = 10.1 kpc (for Rq = 8 kpc), and in this direction the Perseus arm shock (which defines 
where the arm is) is 9.3 kpc from the center. By assuming the arm is a logarithmic spiral 
with pitch angle 12°, one finds the intersection of a 10.1 kpc radius circle with the shock at a 
current position of ~ 5° (here is galactocentric azimuth, defined as zero from the Galactic 
center to the Sun and related to longitude £ by sin0 = (r/R) sin£; hence (p is positive in the 
2'^'^ quadrant of longitude). Thus the stars have since migrated some A0 ~ 30° beyond the 
arm (their formation place) to their current position R, (p = 10.1 kpc, 35°. How long has 
it taken them to move this 30° angular distance? For "flat" circular rotation the angular 
velocity of CTB 102 is = 220/10.1 = 21.8 km s~^ kpc^^. A reasonable pattern speed 
for t he Perseus arm is Qp = 15 km s^^ kpc~^, a mean between mode rn estimates near 20 
(e.g. iMartos et al.ll2004l ) and older estimates near 11 (e.g. lRobertslll972l ) for the Milky Way's 
spiral pattern. Then the time it took the system to migrate A0 from the arm at angular rate 
Q — Qp (the relative angular velocity with respect to the arm) is t = A0/ {Q — Qp) = 75 Myr. 
If the massive stars in CTB 102 formed from the compression of the shock, then by now 
the only ones left would be those that live longer than 75 Myr, or stars of < 7 Mq (B5V 
types, for example). A cluster of hundreds of such stars would be needed to power CTB 102 
[U = 112 pc cm~^; see Sec. 16.11) . It is highly unlikely that CTB 102 contains stars later 
than B5 only. It is more likely that CTB 102 lies near to the arm/shock, which is 4.3 kpc 
distant in this direction. 

The uncertainty in this estimate is ±20%, from the variation of best-fitting models 
of the same direction and in models fitted to several immediately adjacent directions (see 



Foster fc MacWilliamsl l2006l ). This agrees with stellar distances of the nearby optically- 



briUiant H II regions Sh2-124 (£, b = 94?6, -1?5, Vco = -^^ km s-\ r = 3.6 ± 0.6 kpc) 
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and Sh2-132 {i, b = 103°, -0?8, Vco = -49 km s-\ r = 3.5 ± 0.9 kpc), both of which 
are Perseus arm objects (^^CO/H I-based velocities and spect rophotometric dis tances to 100 
Sharpless objects in the second quadrant are forthcoming in Foster et al.ll2009l ). 



We note that CTB 102 and filaments may instead be kinematically associated with the 
Cygnus spiral arm, participating in an extensive star-forrnation group along with H II region 
NRAO 655 ( VrsR 70 km s-^ IPoster fc Routledeel boOll ) and SNR 3C434.1 {Vlsr ^ 



-80 km s 



Foster et al.ll2004j ). This would extend its distance estimate here by a factor 
of 1.6, since the mean stellar distance to Cygnus arm H II regions between 90° < i < 100° 
(Sh2-121, BFS 8, Sh2-128, BPS 10 & DA 568) is r ~ 7.0 kpc. 



6. Radio Continuum Analysis 

Radio continuum measurements of H II regions can provide a lower limit on the total 
number of ionizing photons (A'"^,) and estimates of the average excitation parameter (U). 
Assuming an opticall y thin, spherical, constant-density H II region, the equations used are 



(IRudolph et al. 



19961 ): 



Nl = 7.5 X 10^3^,ciV°-iT-°-^^ s-^ (3) 

U = 1.335y3,i2/3j,l/30y0.116 ^^-2 

where u is the frequency in GHz, 5*;^ is the flux density measured at frequency u in mJy, d 
is the distance to the source in kpc, and Tg is the electron temperature in units of 10^ K. 
Parameter values used are u = 1.42 GHz, d = 4.3 kpc and Te = 10000 K. Since none of 
the parameters are ver y sensitive to Tp, and the uncertainties in Table [3] are large, a general 



value of 10^ K is used ( iRudolph et al.lll996l ). 



The flux density for the considered region is obtained from the CGPS continuum image 
by using the Dominion Radio Astrophysical Observatory (DRAO) software IMVIEW. A 
polygon incorporating the region of interest is drawn by eye. The polygon's perimeter deflnes 
the background, and a twisted plane is fltted to estimate the background. The flux density 
with the estimated background subtracted is then used to compute the parameters Nl and U. 
The biggest contribution to the uncertainty is where to draw the polygon around the region 
of interest. To estimate this uncertainty, each flux measurement is done ten times, each time 
with a new polygon drawn. The average value is used to compute the estimated parameters, 
and the standard deviation is used to estimate the uncertainty in the flux measurement. 
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6.1. CTB 102 Central Region 

A representative polygon around the region of interest is shown in Figure [6l top. The 
flux density is measured ten times with results ranging between 28.9 and 36.9 Jy. After 
subtraction of the flux density from the foreground source WB 43, the average value and la 
uncertainty in the flux density is 31.3 ± 2.5 Jy, which gives a value of = (4.5 ± 1.8) x 
10^^ s~^ and U = 112 ± 15 pc cm~^. Almost all the uncertainty in these valu es comes 



from the uncertainty in the distance. Comparing the values (logA'^i;, = 49.7) to iPanagia 



( I1973I ) sugges ts an ionizing flux a nd excitation parameter corresponding to a 05 V class star. 



Comparing to lVacca et al.l (119961 ) suggests an 04 V star. However, the 30 K contour in Figure 
E] shows the brightest parts of CTB 102 being along an arc stretching from CTB102p towards 
fllament 7 (i.e. towards smaller Galactic longitude and latitude). Given this continuum flux 
distribution, it is more likely that the ionizing flux comes from several stars. In summary, 
the estimates together with the morphology of the 1.42 GHz continuum emission, suggest 
that the CTB 102 complex is powered by at least several late type O stars, maybe even an 
early type O star (since is a lower limit). 



6.2. Filament 10 



Filament 10 looks strikingly circular in the 1.42 GHz continuum image (see Figure O 
bottom) with an angular diameter of ~ 5', yielding a physical size of ~ 6 pc at a distance of 
4.3 kpc. The intensity proflle, with a rapid fall-off at the edge, strongly suggests the region is 
ionization bound. As such this filament would not contribute ionizing photons to the larger 
complex. Drawing a circular polygon around Filament 10 (average background ~ 11.1 
K) yields 0.46 ± 0.03 Jy for the flux density (same method for uncertainty estimate as for 
the whole complex), giving Nl = (6.6 ± 2.6) x 10^^ s "^ (logiV^ = 4 7.8) and U = 28± 4 pc 



cm~^. Comparisons suggest an 09.5 V (|Panagialll973l ) or an BO. 5 V (IVacca et al.lll996l ) star 
ionizing the region. For this fllament, a single star or a compact cluster of early type B stars 
are ionizing the probably homogeneous ISM, given the circular radio continuum proflle. 



7. Discussion 

This study was partly motivated by the appearance of the CTB 102 region in the 1.42 
GHz continuum mosaic (Figured]). Looking at Figure [31 one can make out the region looking 
like a cone-shaped bubble. The radio brightness is concentrated along a ridge following an 
arc with a maximum at CTB102p. The arc goes towards fllaments (in turn) 7, 6, 5, 4 and 
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continues towards greater Galactic latitude, with faint filamentary structures around b = +4° 
eventually bending back towards greater Galactic longitude and then down. Within the arc, 
there appears to be a "cleared out" region extending from the central region of CTB 102 
towards greater Galactic latitude and smaller Galactic longitude. The observed velocities for 
the filaments along this arc all fall within < 11 km s~^ of the velocity of CTB102p, indicating 
that they at least are in the same part of the Galaxy. Taking the velocities and 1.42 GHz 
continuum image together, it is likely that the CTB 102 complex includes filaments 7, 6, 5 
and 4. Filament 2 is probably not part of the region, since its velocity deviates by 16 km s~^, 
and no obvious association with the other filaments is seen in the radio continuum data. As 
will be discussed below, our model shows that filament 2 could not be powered by CTB 102. 

With the inclusion of these filaments, the angular size of the CTB 102 H II region is 
~ 80', which corresponds to a physical size of ~ 100 pc at a distance of 4.3 kpc. Inclusion 
of filaments 8, 9, and 10, all within < 6 km s^^ of CTB102p with apparent connections to 
CTB 102 in the 1.42 GHz image, increases the angular size of the H II to ~ 100' or ~ 130 
pc at 4.3 kpc. The objects KR 4, KR 6 and NRAO 652 have no obvious connection to the 
region in the radio data, and are probably H II regions which happen to be in the same part 
of the Galaxy as CTB 102. 

If the assumption of the filaments being part of the structure is indeed correct, what 
sort of process would create such a big H II region? The typical sound speed in a H II region 
is on the order of 10 k m s~^, which is also t he typical velocity associated with blister regions 



and champagne fiows (jTenorio-Tagldll982i ). This velocity corresponds to ~ 10 pc/(10^ yr). 
The velocity is consistent with the observed velocities in Table [2l However, to reach a radius 
of 50 — 65 pc {2R =~ 100 or ~ 130 pc) would require a time scale of 5 to 6.5 million years 
(assuming constant expansion velocity). This time scale is lon ger than the main-sequence 



life time of a 05 V star, 3.7 x lO'' years (ISchaller et al.lll992h . but comparing CGPS 408 



MHz and 1.42 GHz show no sign of a supernova in the CTB 102 region. Also, CTB 102 does 
not appear like a blister in the 1.42 GHz continuum image. 

Rather, an alternative interpretation is the whole CTB 102 complex with its filaments 
is a combination of an H II region with a bubble/chimney structure. A stellar wind from 
massive a O star(s) creates a bubble in the interstellar medium (ISM). Towards a region 
of the Galaxy with lower density (higher Galactic latitude) the expanding bubble has its 
top "blown off", creating the possible chimney. Towards regions with higher densities, the 
expanding stellar wind bubble sweep up a shell of the ISM. The shell and left over neutral 
material is exposed to ultraviolet radiation from the O star(s), gets ionized and forms the 
H II region. 



What would be expected if this interpretation is correct? The H I structure would show 
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a cavity in H I. In this picture, the expanding bubble sweeps up a shell of neutral material. 
This shell would be ionized from the inside of the bubble and if viewed through the edge, 
the increased path length would cause the shell to appear as neutral and ionized filaments. 
These H I and H II filaments would have observed velocities which do not deviate too much 
from each other and from the velocity of the central region. Also, the ionized filaments would 
be expected to be located on the "inside" of the neutral filaments. 

Support for the idea comes from viewing the H I structure in the region. Figure [7] shows 
the averaged velocity slices (of width 2.5 km s~^) of the CGPS H I line data, overlayed with 
contours (15 and 30 K) indicating where the 1.42 GHz continuum flux is located in relation 
to the H I. A cleared out region is most easily seen in the final three velocity frames (—60.8 
to —65.8 km s~^), extending from the central region of CTB 102 towards the Galactic north. 
This cavity in H I disappears beyond ~ —72 km s~^. The Galaxy contains gas out to ~ —100 
km s~^ in this direction (see Figure HJ center). The H I emission is mostly concentrated along 
an arc extending from north of CTB 102, wrapping around the "edges" of the H II region 
and extending more faintly up next to filaments 7, 6 and 4. These H I filaments (seen in 
the final three velocity frames, —60.8 to —65.8 km s~^) are located on the "outside" of the 
corresponding H II filaments and they correspond roughly in velocity to the RRL results in 
Table [21 In the first velocity frame, —53.4 km s~^, a filament of H I is spatially coincident 
with filament 5 (observed Vlsr = —51.9 km s~^). 

Figure [S] shows a MSX A-band image (8.3 fim) towards CTB 102. The contour corre- 
sponds to an averaged CGPS H I brightness temperature level of 70 K from the —55.9 km 
s~^ velocity frame in Figure [3 The MSX A-band traces PAH emission from neutral ma- 
terial being bombarded with ultraviolet photons. Figure [8] clearly shows the morphological 
association between structures seen in the H I arc discussed above, neutral material in a UV 
radiation field and the 1.42 GHz continuum emission. This is best demonstrated by the 8.3 
lim emission tracing a thin filamentary H I structure which extends from the main H I arc 
towards CTB102p (i.e. towards greater Galactic longitude and latitude). This thin filament 
runs parallel to the bright arc seen in the 1.42 GHz continuum image (traced by the 30 K 
contour in Figure [6]), something which would be expected if the interpretation of a neutral 
shell being ionized from slightly above the bright arc is correct. 

Can ionizing photons reach the filaments? To reach the inside of the proposed swept 
up shell at filament 6, ionizing photons have to travel an angular distance of ~ 50' from 
CTB102p, corresponding to ~ 60 pc at 4.3 kpc. In the bubble scenario, the inside of the 
bubble is filled with stellar wind, a hot, low density plasma. Absorption in such a plasma is 
negligible, and the purely geometrically diluted radiation field could easily maintain a high 
ionization fraction in the ionized inside of the swept up shell. The same is true for the other 
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filaments we consider part of tlie sliell structure. For filament 2, at a distance of ~ 180 
pc from CTB102p (~ 140'), a purely geometrically diluted radiation field could maintain a 
high ionization fraction if the stellar wind bubble extended all the way out to filament 2. 
The proposed bubble does not extend that far (see below), and inserting any realistic ISM 
between the proposed bubble and filament 2 causes significant absorption in the ISM, thus 
leaving no photons from CTB 102 to ionize filament 2. 

A simpl e analytical model of the develo pment of a s tellar wind bubble was given by 



Castor et al.l (Il975l ). Using the formulation in iKwokl (120071 ). the radius (R), shell expansion 



velocity [R) and shell fractional thickness [AR/R) of such a region are given by: 

(t/10^ yr)^ pc (5) 

(t/10^ yr)"^ km (6) 

(7) 

where E is the energy loss rate from the star, /i is the mean molecular weight of the medium 
the shell is expanding into (/i = 1.4 for neutral ISM, /i = 0.62 for a fully ionized ISM), no is 
the ISM density and is the sound speed inside the bubble. 

We estimate E using a typical mass loss rate of M = 10^^ Mq/jt and a stellar wind 
velocity of V = 2000 km s~^ The so und speed is an = 10 km s~^ and the ISM density is 
uq = 1 cm~^ (mean ISM density, from Basu et al. 19991 ). If /i = 1.4 it would take 3.8 x 10^ 



R = 


28 


( E/W^ erg s"^ 


I /urio/l cm^'^ 


R = 


16 


/ E/IO^^ erg s"^ 


I fmo/1 cm~3 


AR 






R 




yRj 



years to reach the position of filament 6 from CTB102p (60 pc), by that time the expansion 
velocity of the shell will be 9 km s~^. The other extreme, if /i = 0.62, a radius of 60 
pc would be reached in 2.8 x 10^ years, with a expansion velocity at that time of 12 km 
s^^. The expansion velocity for both scenarios is fairly insensitive to large t, and reaches 
order of magnitude ~ 10 km s~^ after roughly 2.5 million years. The estimated thickness 
of the proposed shell from the 1.42 GHz continuum image and the H I line data is ~ 6', 
corresponding to a fractional thickness of ~ 12%. The time it takes the model shell to reach 
such a fractional thickness is 0.9 x 10® years (if fi = 1.4) or 1.3 x 10® years (if fi = 0.62). 
If instead no = 10 cm~^ (the rest of the parameters remain the same), time to reach 60 
pc increases by a factor of ~ 4. This suggests that if this mechanism is responsible for the 
observed filaments, the ISM in the direction of decreasing Galactic latitude cannot be more 
than a few particles per cubic centimeter, or the timescale exceeds the main-sequence life 
time of the probable star(s). 
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In the basic model for the expanding swept-up shell, the shell density would be (iKwok 



20071 ): 



Ps= {—] Po cm 3 (8) 

Estimates of the column density test this picture where the observed filaments are due to 
our line of sight through the edge of the shell. The maximum path length s looking through 
the edge of a uniform density spherical shell of radius R and thickness AR would be: 



s = 2R\J2 (AR/R) + {AR/RY (9) 

Estimating the column density towards the filaments is done by integrating the CGPS H I 
line data. The single channel (width 0.82 km s~^) rms noise is 3.3 K, corresponding to a 
single channel rms noise in the column density of 4.9 x 10^^ cm~^. For filament 6, i? = 60 pc 
and AR/R = 0.12, which gives s ~ 60 pc. This is on the order of the radius of the proposed 
bubble in this direction, so this is probably an upper limit in path length through the edge of 
any shell associated with CTB 102. The column density through filament 6 is ~ 10^^ cm^^ 
(making the rms noise negligible). The uniform density needed in a shell to produce such a 
column density is ~ 5 cm~^. In the model, the time required to reach this shell density (with 
uq = 1 cm~^) is 0.6 X 10^ years. Given the large uncertainties of the parameters involved, in 
R (from the distance uncertainty), AR (uncertainty in telling where the filaments start and 
end) and column density (mainly determining the background), the uncertainty in the shell 
density can be as large as ~ 60%. Within such a large uncertainty, the model can account 
for observed column densities within a time frame set by the main-sequence life time of a 
massive star. However, it cannot be ruled out that the observed filaments are higher density 
regions in the ISM that happen to be exposed to ionizing photons from the star(s) powering 
the H II region. Even if that is the case, the medium between the filaments and the powering 
source would most likely be low density material with a high ionization fraction to allow for 
enough photons to reach the filaments. Such a medium is consistent with the H I line data 
and the 1.42 GHz continuum image. 



The ICastor et al.l (119751 ) model is very basic and does not take into account multiple 
stellar winds nor non-uniform ISM, but the estimates indicate that the bubble model is 
a plausible explanation of the size of CTB 102, allowing for observed size, velocities and 
fractional shell thickness in a time scale of 1 — 2 x 10^ years. That is less than the main- 
sequence life time of the star(s) that is probably powering the region. For the same reason, 
it is unlikely that the proposed bubble would reach all the way out to filament 2. The 
overall morphology of CTB 102 indicates that the proposed H II region/bubble /chimney is 
a major feature in the Perseus Galactic arm, powerful enough to disrupt the ISM and clear 
out a ~ 100 — 130 pc region. The appearance of CTB 102 in the H I with a cleared out 
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regio n makes the CTB 102 com plex comparable in size to the W4 supe rbubble in the Pers eus 
arm (INormandeau et al.lll996l ). The W4 superbubble is modeled by iBasu et al.l (119991 ) as 
expanding in a stratified atmosphere. An association of nine 0-type stars (IC 1805) is 
probably the reason for the cavity in H I and the ionization structure of the W4 H II region. 
CTB 102 shares similarities with W4 in large scale H I distribution and ionization structure. 



8. Summary 

We have obtained new RRL observations of CTB 102 that show that the filamentary 
structure surrounding the central region is physically associated with the central region. We 
provide the first ever distance estimate for this H II region, 4.3 kpc. We argue that the best 
explanation for the size and appearance of the whole complex is that it is a large H II region 
combined with a wind-blown interstellar bubble/chimeny structure. MSX, H I and 1.42 GHz 
observational data are consistent with this view. 
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Fig. 1. — CGPS 1.42 GHz continuum image with GET pointings (see Tabled]) labelled and 
displayed as circles. The size of the circles correspond to the GET beam size, 2.5'. The 
contours correspond to 15 and 30 K brightness temperature levels. 



-17- 




-100 -50 50 
Vlsr (km s-') 



Fig. 2. — Reduced spectra (solid) and Gaussian fits (dashed) for eacli observation, plotted 
on the same vertical and horizontal scale. The number in parenthesis is the scaling factor by 
which the spectrum was multiplied in order to fit the scale. The dash-dotted lines correspond 
to the 3(7 and 5a noise levels. The vertical dotted line is the central velocity from the Gaussian 
fit. 
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Fig. 3. — As Figure [T] but showing the observed velocities for each position. 
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Fig. 4. — Three views of the large-scale H I environment near CTB 102, taken from the 36' 
resolution 26-m DRAO H I survey. (Top) {b, V) (latitude-velocity) centered on i = 93?6, 
(Center) Tb {i, V) (longitude-velocity) centered onb = +2?8, (Bottom) Tb {i, b) (longitude- 
latitude) centered at Vlsr = —58 km s~^. The position i, b, V of CTB 102 is marked with 
the star symbol in each slice. Brightness-temperature contours in each panel are at 40, 54, 
65 and 80 K. CTB 102 appears associated with H 1 that could either be a high- latitude 
extension of the Perseus spiral arm (centered at —44 km s~^) pushed to negative velocities 
by the "rolling" motion in the arm, or part of a low- velocity extension of H I from the Cygnus 
spiral arm (—68 km s~^), either of which could be part of a very large H 1 bubble expanding 
away from its parent arm. 
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Fig. 5. — Model velocity field t owards I = 92?4, b = — 0?2, from the best-fitting model 



of H I density and kinematics of iFoster &: MacWilliamd (120061 ) . The circular velocity field 
only (no streaming motions from density waves) is shown by the thick line; with streaming 
added, the result is the thin line which includes the shock at 4.2 kpc. Based on CTB 10 2's 
velocity (corrected for the "rolling" gradient in the Perseus arm, as in [Foster et al.ll2009l ) of 
^ls'r ~ ~52 km s~^ is found at either 4.3 kpc or 5.8 kpc along the line-of-sight (assuming it 
is a Perseus arm object; see text). 
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Fig. 6. — Two CGPS 1.42 GHz continuum images of CTB 102 overlayed by contours cor- 
responding to 15 and 30 K. (Top) The dashed polygon is a representative of the polygons 
used in Section 16. 1[ (Bottom) The circle indicates the area used for estimates for filament 
10 in Section [621 
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Fig. 7. — CGPS H I line image, averaged over three velocity channels to a width of 2.5 km s~^, 
centered on the velocity shown in the frame. The white contours correspond to 1.42 GHz 
continuum brightness temperature levels of 15 and 30 K. The white circles indicate the 
positions of the GET pointings for filaments 4, 5, 6 and 7 (from Figure [1]). This figure 
is reduced due to size constraints. For the full version, see the published article in The 
Astrophysical Journal. 
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Fig. 8. — MSX A-band image (8.3 /im). The contour corresponds to an averaged CGPS H I 
brightness temperature level of 70 K from the —55.9 km s~^ velocity frame in Figure [71 
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Table 1. Observations. 



Filament / 


RA 


Dec 




h 


Integration time 


T-^ system 


Object 


h m s 


o / // 


O 


O 


s 


K 


CTB102p 


21 12 28.9 


52 32 23 


93.115 


2.835 


3600 


23.7 


2 


21 01 32.5 


50 51 17 


90.725 


2.955 


30000 


21.9 


4 


21 05 19.1 


52 17 01 


92.185 


3.465 


14400 


21.3 


5 


21 05 56.4 


51 52 53 


91.950 


3.125 


15000 


21.0 


6 


21 07 14.6 


52 07 03 


92.260 


3.135 


10200 


21.1 


7 


21 10 48.6 


51 55 51 


92.495 


2.605 


3000 


21.9 


8 


21 15 39.9 


52 22 03 


93.325 


2.365 


7200 


22.4 


9 


21 13 59.6 


52 15 09 


93.065 


2.470 


2400 


23.0 


10 


21 13 49.2 


52 27 25 


93.195 


2.630 


600 


23.5 


KR4 


21 16 24.3 


52 48 02 


93.715 


2.585 


3000 


22.1 


KR6 


21 19 01.0 


53 22 26 


94.400 


2.705 


6600 


21.5 


NRAO 652 


21 15 20.1 


51 52 56 


92.940 


2.065 


6600 


23.5 


WB 43 


21 09 21.6 


52 22 22 


92.668 


3.069 


600 


23.9 



Note. — Positions are J2000. Integration times and average system temperatures 
are for the GET observations. 
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Table 2. Observed line parameters. 



Filament / 


Ti 


AT 

'—^-^ rms 






1^ 




Object 


(mK) 


(mK) 


(km s~^) 


(km s~-^) 


(k, 




CTB102p 


102.6 ±0.4 


1.2 


-62.66 ±0.05 


18.34 ±0.08 






2 


1.7±0.1 


0.3 


-78.4 ±0.6 


13.9 ± 1.1 




16 


4 


2.6 ±0.1 


0.6 


-67.2 ±1.3 


30.6 ±2.2 




5 


5 


3.5 ±0.2 


0.7 


-51.9 ±0.9 


21.0± 1.5 




11 


6 


3.5 ±0.2 


0.6 


-67.0 ±0.8 


25.4 ± 1.4 




4 


7 


6.7 ±0.3 


1.2 


-57.7 ±0.8 


22.8 ± 1.4 




5 


8 


4.4 ±0.2 


0.7 


-56.4 ±0.7 


21.2 ±1.2 




6 


9 


10.7 ±0.3 


1.1 


-66.3 ±0.5 


24.8 ±0.8 




4 


10 


29.0 ±0.5 


1.5 


-65.0 ±0.2 


19.7 ±0.4 




2 


KR4 


12.1 ±0.3 


0.9 


-67.8 ±0.4 


24.5 ± 0.6 




5 


KR6 


4.7 ±0.2 


0.7 


-56.1 ±0.5 


16.0 ±0.9 




7 


NRAO 652 


5.7±0.2 


0.7 


-68.7 ±0.5 


26.0 ±0.9 




6 


WB 43 


68.9 ±0.6 


1.7 


-0.8 ±0.1 


15.3 ±0.2 




62 



Note. — Spectral parameters for each filament/object observed. Line amplitude 
(T;), central velocity (Vlsr) and FWHM (AV) are obtained by a Gaussian fit to 
the radio recombination line after baseline subtraction, regridding and smoothing. 
The noise, ATj-ms, is obtained by considering regions on both sides of the line. 
Both T; and ATrms are given in brightness temperature units. The uncertainties 
correspond to la. The \ V — Vref\ column displays the absolute difference in Vlsr 
for each observed filament/object, compared to CTB102p and rounded to the 
nearest integer. 
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Table 3. Derived parameters. 



Filament / 


Td 


Te 


Object 


(K) 


(K) 


CTB102p 


7400 ± 100 


5200±-° 


2 


4200 ± 600 


3400+^2200 


4 


21000 ± 2900 


ciQnn+3600 


5 


9800 ± 1300 


6000+3^°° 

DUUU_43QQ 


6 


14400 ± 1500 


r,nnn+-''2oo 

dUUU_3gQg 


7 


11600 ± 1300 


5000t^«°° 


8 


10000 ± 1100 


83001^900 


9 


13800 ± 800 


790011?°° 


10 


8600 ± 300 


7700t^S° 


KR4 


13400 ± 600 




KR6 


5600 ± 600 


5800+3^°° 


NRAO 652 


15100 ±1000 




WB 43 


5200 ± 100 


5500t|°° 



Note. — Derived parameters for each fil- 
ament observed. The uncertainties in Td 
are la uncertainties. The listed uncertain- 
ties in Te comes from considering uncer- 
tainty in the background estimate, 0.5 K. 
The total uncertainty in Te are by far dom- 
inated by this background uncertainty, the 
other two contributions (from T; and AV) 
are both ~ 10% of the dominant back- 
ground contribution. 



